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Abstract 

It is nontrivial whether the average size of a ring polymer should become smaller or 
larger under a topological constraint. Making use of some knot invariants, we evaluate 
numerically the mean square radius of gyration for ring polymers having a fixed knot 
type, where the ring polymers are given by self-avoiding polygons consisting of freely- 
jointed hard cylinders. We obtain plots of the gyration radius versus the number of 
polygonal nodes for the trivial, trefoil and figure-eight knots. We discuss possible 
asymptotic behaviors of the gyration radius under the topological constraint. In the 
asymptotic limit, the size of a ring polymer with a given knot is larger than that 
of no topological constraint when the polymer is thin, and the effective expansion 
becomes weak when the polymer is thick enough. 
PACS number(s): 36.20.-r, 61.41. +e, 05.40.Fb 
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1. Introduction 



The effect of a topological constraint should be nontrivial on physical quantities of a ring 
polymer such as the size of the ring polymer. Once a ring polymer is formed, its topological 
state, which is given by a knot, is fixed. However, it has not been established how to for- 
mulate the topological constraint on the ring polymer in terms of analytic methods. On the 
other hand, several numerical simulations have been performed, investigating some statisti- 
cal properties of ring polymers under topological constraints^^. Through the simulations, 
it has been found that a topological constraint may severely restrict the available degrees of 
freedom in the configuration space of a ring polymer, and can be significant in its physical 
properties. 

Recently, DNA knots are synthesized in experiments, and then separated into various 
knot types by the agarose-gel-electrophoresis techniqueBHll. Under the electric field, the 
charged macromolecules move through the network of the gel, and the migration rates should 
depend on their sizes, shapes and charges. It is remarkable that the electrophoretic mobility 
of a circular DNA depends also on its knot type. It is observed that the more complicated 
the DNA knot is, the higher its mobility. The fact could be related to a common belief that 
the gyration radius of a knotted DNA should depend on its knot type. 

In this paper, we discuss the mean square radius of gyration of circular polymers having 
a fixed knot type. We consider the question how the size of a ring polymer should depend on 
the length N under the topological constraint. Here the length N corresponds to the number 
of polygonal nodes, when we model the ring polymers by some self-avoiding polygons. We 
find that the effect of the topological constraint is not trivial. In fact, the mean square radius 
of gyration for ring polymers under the topological constraint can be smaller or larger than 
that of the ring polymers under no topological constraint, as we shall see later. We study the 
question through numerical simulations, making use of some knot invariants. We construct 
an ensemble of ring polymers consisting of freely jointed hard cylinders, and then evaluate 
the mean square radius of the gyration for the ring polymers that have the given knot type. 
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2. Cylindrical self-avoiding polygons 




Let us explain the model of ring polymers consisting of freely jointed hard cylindersoia. 
The segments are given by hard cylinders with radius r. They are "hard" in the sense that 
there is no overlap allowed between any pair of non-adjacent cylindrical segments, while 
adjacent segments can overlap: there is no constraint on any pair of adjacent cylinders. The 
model was first introduced in Ref.i with the Monte-Carlo algorithm using the "hedgehog" 
configurations of polygons. Recently, another method was introduced for constructing the 
cylindrical self-avoiding polygons^!. It is based on the algorithm of ring-dimerizationi for 
the rod-bead model of self-avoiding polygons. We call the new algorithm the cylindrical 
ring-dimerization method, or the dimerization method, for short. All the cylindrical self- 
avoiding polygons in this paper are constructed by the dimerizaion method. We note that 
the algorithm of chain dimerization is quite effective for constructing long self-avoiding walks 
on off-latticeEl. 

The most useful property of the model of cylindrical self-avoiding polygons is that it 
has the two parameters: the cylinder radius r and the number N of segments. The two 
parameters are important in the theoretical explanation of the experimental data of DNA 
knotsi: the cylinder radius corresponds to the effective radius of negatively charged DNAs 
surrounded by the clouds of counter ions. Thus, the model is closely related to the wormlike 
chain model for polyelectrolytes in electrolyte solutions. Furthermore, it has been shown 
that the cylinder radius r plays an important role in the random knotting probability^. 

3. Gyration radius of cylindrical self-avoiding polygons 

Now let us consider the mean square radius of gyration R 2 for the cylindrical self-avoiding 
polygons. It is defined by 
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Here R n is the position vector of the n-th monomer and iV is the number of nodes in the 
cylindrical self-avoiding polygon. In numerical simulations, the statistical average < • > 
is given by the average over M samples of polygons. We take M = 10 4 in the paper: 
we construct 10 4 self-avoiding polygons of N segments of cylinders with radius r by the 
dimerization method. 

We have calculated numerical estimates of the the mean square radius of gyration R 2 for 
the cylindrical self-avoiding polygons with the different values of radius r for several different 
numbers N of polygonal nodes upto 1000. We consider eighteen numbers from 20 to 1000 
for the number N of polygonal nodes, and fifteen different values from 0.0 to 0.07 for the 
cylinder radius r. Here we note that the gyration radii for the rod-bead polygons have been 
estimated in RefsBS. 

4. Method for selecting polygons with the same given knot 

Let us describe the method for selecting such polygons that have a given knot type in our 
simulation. For a given knot K, we enumerate the number Mr of such polygons out of the 
M polygons that have the same set of values of some knot invariants for the knot type K. 
We employ two knot invariants, the determinant A#(— 1) of knot and the Vassiliev-type 
invariant v 2 (K) of the second degree, as the tool for detecting the knot type of a given 
polygoniiii. 

The number depends not only on the knot type but also on the step number N and 
cylinder radius r. Let us recall the knotting probability for cylindrical self-avoiding polygons 
consisting of cylinders with the radius r0. We denote by P tr i V (N,r) the probability of a 
cylindrical self-avoiding polygon of N nodes with cylinder radius r being a trivial knot. 
Then, it is given by 

Ptnv(N,r) = exp(-N/N c (r)). (2) 
Here N c (r) is called the characteristic length of random knotting and can be approximated 
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by an exponential function of r: N c (r) = N c (0) exp(7r), where 7 is a constantl 



5. Gyration radius of ring polymers under a topological constraint 

Let us discuss the mean square radius of gyration for cylindrical self-avoiding polygons with 
a fixed knot type. We evaluate the expected value of the gyration radius for a given knot 
K by the mean value of the square radii of gyration for the Mk polygons with the knot K. 
Then, the mean square radius of gyration R 2 K for the knot type K is given by 

1 Mx 

R K = Tl- R K,ii (3) 
M K i=1 

where R 2 Ki denotes the gyration radius of the i-th cylindrical self-avoiding polygon that 
has the knot type K, in the set of M polygons. For the rod-bead model of self-avoiding 
polygons, R\ has been evaluated for some different knots0. 

Let us discuss our data of numerical estimates of R 2 K . For trivial, trefoil and figure-eight 
knots, we find that the mean square radius of gyration increases monotonically with respect 
to the number N of polygonal nodes. However, the ratio R K /R 2 is not constant with respect 
to N. The graphs of the ratio R 2 K /R 2 against the step number N are plotted in Fig. 1 for 
trivial and trefoil knots. Here the cylinder radius is given by 0.005, i.e., the diameter is 
given by 0.01. Here we recall that R 2 denotes the mean square radius of gyration for the 
self- avoiding polygons with all possible knot types. In terms of R K , R 2 is given by the 
following: R 2 = J2k m kR k / M - We reca11 that the number M K gives different values for 
different step numbers N or different knots K. When the knot K is complicated, Mk can 
be very small and it may give a poor statistics to R\. 

Let us discuss the plots of the trivial knot shown in Fig. 1. The ratio R 2 triv jR 2 at N = 21 
is almost given by 1.0. This is consistent with the fact that trivial knots are dominant when 
iV is small. Here we recall that the probability of being a trivial knot is given by eq. (0) . 
On the other hand, the ratio R 2 riv /R 2 increases with respect to the step number N. Thus, 
the size of the ring polymer enlarges under the topological constraint of being a trivial knot. 
The topological constraint gives an effective swelling effect, in this case. 
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Let us consider the case of trefoil knot. The ratio R 2 tre jR 2 is not always larger than 
1.0. In fact, it is smaller than 1.0 when iV < 200. The ratio R 2 re /R 2 is given by about 0.7 
when N = 21. Thus, the topological constraint gives an effective shrinking effect on the ring 
polymer. When N > 300 or 400, the ratio R 2 re /R 2 becomes larger than 1.0 for r = 0.005. 
Thus, the topological constraint makes the ring polymer enlarge for large N. However, when 
r > 0.03, the ratio R 2 re /R 2 becomes smaller than 1.0 even for N = 1000, as shown in Fig. 
2. 

In Figs. 1 and 2, we see that the difference among the ratios of the gyration radii for 
trivial, trefoil and figure-eight knots is much more clear in the small- N region than in the 
large- N region. The dependence of the gyration radius of the ring polymer on its knot type 
could be more significant for small N than for large N. The effective shrinking of a thick 
and finite ring-polymer with a nontrivial knot might be associated with the concept of ideal 
knots or tight knotsHS. 

The fitting curves in Figs. 1 and 2 are given by 

?k= lK (l-6 K eM-VKN)) . (4) 

Here, the constants jk, &k and tjk are fitting parameters. We see that the fitting curves 
in Figs. 1 and 2 are consistent with all the numerical estimates of R\ in the range from 
N = 21 to 1000. Thus, the formula (|j) effectively describes the finite-size behaviour of R 2 K , 
although there is no apriori reason for assuming it. For instance, the formula (f|) can be 
not appropriate as an asymptotic expansion. In §6, we shall introduce another formula to 
discuss the possible asymptotic behaviors of R\. 

6. Asymptotic behaviors of R 2 K 

Let us discuss possible asymptotic behaviors of the gyration radius of the ring polymer 
under a topological constraint. We may assume that when N is very large, R 2 K can be 
approximated by 
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R\ = A K N 2UK (l + B k N- Ak + 0(l/N)) . (5) 

The expansion is consistent with renormalization group arguments, and hence it should be 
valid for the case of asymptotically large N. The exponent Vk and the amplitude Ak can 
be evaluated by applying the formula (|5J) to the numerical data of R 2 K for large values of N. 
Here we note that the expansion is not effective for small N. In fact, when N < 200, it 
does not give any good fitting curves to the data in Figs. 1 and 2. 

Let us now discuss the exponent Vk- We have analyzed the plots of R 2 K /R 2 given in 
Fig. 1 for trivial and trefoil knots by the least-sqaure method with respect to the following 
formula: R 2 K /R 2 = {A K jA)N VK ~ v (l + (B K - B)N~ A + 0{l/N)\. Here we have assumed 
Ax = A = 0.5, which is consistent with the scaling expansion. Then, we obtain the following 
estimates: for trivial knot, v K -v = 0.005±0.103, A K /A = 1.18±0.99, B K -B = -1.7±4.0; 
for trefoil knot, u K -u = 0.009 ±0.090, A K /A = 1.18 ±0.85, B K -B = -3.5 ±3.0. To each 
of the knots, we have applied the formula (|5|) to the eight points with iV > 300 in Fig. 1. 
The x 2 values are given by 4.3 and 6.2 for trivial and trefoil knots, respectively. For other 
values of radius r, we have similar results. Thus, for trivial and trefoil knots, the exponent 
Vk should agree with the exponent v of the mean square radius of gyration R 2 , within the 
error bars. 

There are also other evidences supporting vk = v, i-e., Vk = 0.588, for trivial and trefoil 
knots. In fact, the plots of the ratio R 2 K /R 2 versus the number N in Figs. 1 and 2 are 
likely to approach some horizontal lines at some large N. It is also the case with some other 
values of cylinder radius r. It is clear particularly for trivial knot. Thus, at least for trivial 
knot, we can easily conclude that the exponent vk should coincide with the exponent v. We 
note that even the fitting formula (^) is consistent with the coincidence of the exponents: 
Vk = v. Thus, all the numerical results obtained in the paper suggest vk = v for trivial 
and trefoil knots. It is also consistent with the observation for the self-avoiding polygon on 
a lattice in Refs.00 that the exponent vk should be independent of the knot type. 

Let us now consider the amplitude A K of the asymptotic expansion (H). In Fig. 
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3, the numerical plots of the ratio Ak/A versus the radius r are shown for triv- 
ial and trefoil knots. The ratio is evaluated by the following formula: R 2 K /R 2 = 
(Ak/A) (l + (B K - B)N~ A + 0(l/iV)) Here, we have assumed v K = v and A K = A=0.5. 
Furtheremore, we have applied it only to the data with iV > 300. 

When r is small, the ratio Ak/A in Fig. 3 is larger than 1.0 for both trivial and trefoil 
knots. It is remarkable that the asymptotic ratio Ak/A can be larger than 1.0 when the ring 
polymer is thin enough. We can easily confirm it also in Fig. 1, observing the increasing 
behavior of the plots of R 2 triv / R 2 versus N. Thus, the topological constraint gives an effective 
expansion to the ring polymer with small radius r. 

Interestingly, we see in Fig. 3 that the ratio Ak/A decreases monotonically with respect 
to the cylinder radius r. One might expect that the ring polymer with larger excluded- 
volume should become much larger. In reality, however, for the ring polymer under the 
topological constraint, the ratio Ak/A decreases if the excluded- volume parameter becomes 
large. It is as if the excluded volume effect could make weaken the effective expansion derived 
from the topological constraint. We notice in Fig. 3 that Ak/A might become close to the 
value 1.0 when r is thick enough. Due to the poor statistics, we can not determine whether 
it really does or not. However, if it does, then it is consistent with the interpretation on the 
lattice model of Refs.0'0 that Ak should be independent of knot type. 

It is quite nontrivial that Ak/A which is valid in the asymptotic expansion can be larger 
than 1.0, and the ratio decreases with respect to the radius r. Let us give one possible 
explanation for it. First, we note that R 2 is given by the average of R\ over all possible 
knots: R 2 = J2k R 2 kPk(N). The fact: A K /A > 1.0 suggests that there are large number 
of knots smaller than the knot K. Second, we recall the finite- N behavior of R 2 re : In Figs. 
1 and 2 we see that the ratio R 2 re /R 2 is much smaller than 1.0 when N is small, while it 
increases with respect to N and finally becomes constant. We consider that when N is small, 
the size of trefoil knot is small due to a finite-size effect, while when N is asymptotically large 
then R\ re /R 2 becomes almost constant. For a given knot K\, if we take N large enough, 
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then the ratio R 2 Ki /R 2 becomes constant with respect to N, while the majority of knots 
possible in iV-noded polygons should be much more complex than the knot K 1 and their 
sizes should be much smaller than that of the knot K\. Thus, R 2 can be smaller than R\ x 
and the ratio A Kl /A can be larger than the value 1.0. 

We remark that R 2 riv /R 2 is always larger than 1.0 both for finite N and asymptotically 
large N. For the finite- N case, it is clear from Figs. 1 and 2 that R 2 riv /R 2 > 1. For the 
asmptotic case, it is suggested from Fig. 3 that A triv /A should be larger than 1.0 for some 
small values of radius r. Thus, the property: R 2 triv /R 2 > 1 should persist in the asymptotic 
limit, as far as the data analysis is concerned. 

Summarizing the numerical results on the asymptotic behaviors, we conclude that the 
topological constraint on a ring polymer gives an effective expansion when the radius r is 
small and also that the expanding effect becomes weaker when the radius r becomes larger, 
and it may vanish when the radius r is large enough. The results suggest that there should 
be a new 'phase transition', where it is controlled by the excluded- volume parameter whether 
the topological constraint gives an effective expansion to the ring polymer or not. 
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FIGURES 

Fig. 1. The ratio R 2 K /R 2 versus the number N of polygonal nodes of the cylindrical self-avoiding 
polygons. Numerical estimates of R 2 riv /R? for r = 0.005 are shown by black circles and those of 
R 2 re /R 2 for r = 0.005 by black triangles. In the inset, the enlarged figure shows the numerical 
estimates of R 2 K /R 2 from N =20 to 100 for the cases of trivial, trefoil and figure-eight knots. 

Fig. 2. The ratio R 2 K /R 2 versus the number N of polygonal nodes of the cylindrical self-avoiding 
polygons. Numerical estimates of R 2 riv /R 2 for r = 0.04 are shown by black circles and those of 
R 2 re /R 2 for r = 0.04 by black triangles. The fitting parameters are given by the follwing: for 
trivial knot, v K -v = -0.0004 ± 0.068, A K /A = 1.04 ±0.59, B K - B = -0.3 ±3.1, and x 2 = 1-1; 
for trefoil knot, v K -v = 0.006 ± 0.101, A K /A = 1.01 ± 0.84, B K -B = -2.0 ± 4.0, and X 2 = 3.2. 
In the inset, the enlarged figure shows the numerical estimates of R 2 K /R 2 from N =20 to 100 for 
the cases of trivial, trefoil and figure-eight knots. 

Fig. 3. The ratio Ak/A versus the cylinder radius r for the cylindrical self-avoiding polygons. 
The values of Ak/A for the trivial knot and the trefoil knot are shown by black circles and black 
triangles, respectively. Each of the black triangles are slightly shifted rightward by one ten-th of 
the value of radius r, for graphical convenience. 
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N : Number of nodes 




N : Number of nodes 




r: cylinder radius 



